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Description 

Estimation Method And System For 
Wpannfa 1 Securities Trading 

The United States Government has certain rights 
to this invention pursuant to grants CCR-91-14042 and 
IRI-92- 12597 awarded by the National Science 
Foundation, and to grant AFOSR-91-0347 awarded by the 
5 U.S. Air Force. 

Background of the Invention 

The invention relates to financial securities 
trading such as, e.g., trading in stocks, bonds and 

0 f inan cial derivative instruments, including futures, 
options and collateralized mortgage obligations. 

In financial securities trading, which includes 
the initial offer for sale, the value of a security may 
be estimated, e.g., based on expected future cash flow. 

5 Such cash flow may depend on variable interest rates, 
for example, and these and other relevant variables may 
be viewed as stochastic variables. 

It is well known that the value of a financial 
security which depends on stochastic variables can be 

0 estimated in terms of a multi- dimensional integral. 
The dimension of such an integral may be very high. 

In collateralized mortgage obligations (CMO) , for 
example, instruments or securities variously called 
tranches, shares, participations, classes or contracts 

5 have cash flows which are determined by dividing and 
distributing the cash flow of an underlying collection 
or pool of mortgages on a monthly basis according to 
pre -specified rules. The present value of a tranche 
can be estimated on the basis of the expected monthly 

0 cash flows over the remaining term of the tranche, and 
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an estimate of the present value of a tranche can be 
r presented as a multi- dimensional integral whose 
dimension is the number of payment periods of the 
tranche. For a typical instrument with a 30 -year term 
5 and with monthly payments, this dimension is 360. 

Usually, such a high- dimensional integral can be 
evaluated only approximately, by numerical integration. 
This involves the generation of points in the domain of 
integration, evaluating or "sampling" the integrand at 

10 the generated points, and combining the resulting 

integrand values, e.g., by averaging. Well known for 
numerical integration in securities trading is the so- 
called Monte Carlo method in which points in the domain 
of integration are generated at random. 

15 With integrands arising in financial securities 

trading, the computational work in combining the 
sampled values is negligible as compared with producing 
the integrand values. Thus, numerical integration 
methods in securities trading may be compared based on 

20 the number of samples required for obtaining a 

sufficiently accurate approximation to the integral. 

Summary Of The Invention 

A preferred method for estimating the value of 

25 a financial security involves numerical integration 

unlike Monte Carlo integration in that an integrand is 
sampled at deterministic points having a low- 
discrepancy property. As compared with the Monte Carlo 
method, significant advantages are realized with 

30 respect to speed, accuracy, and dependability. 

Brief Descriptio n Of The Drawing 

Fig. 1 is a schematic of a programmed computer 
system in accordance with a preferred embodiment of the 
35 invention. 
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Fig. 2 is a graphic representation of performance 
data obtained in computer trial runs with an exemplary 
embodiment of the invention as compared with two Monte 
Carlo computations. 
5 Further included is an Appendix with two computer 

algorithms in "C" source language, respectively for 
computing Sobol points and Hal ton points. For a 
description of C, see B.W. Kernighan et al.. The 
Programming Language C . Prentice-Hall, 1978. 

10 

Detailed Description Of Prefer red Embodiments 

Fig. 1 shows a stored -program computer n con- 
nected to input means 12, e.g., a keyboard, for enter- 
ing financial securities data, and connected to output 

15 means 13, e.g., a visual display device, for displaying 
an estimated value of the financial security. The 
computer 11 includes a working memory M, a low- 
discrepancy deterministic point generator P, an inte- 
grand evaluator E, and an integrand -value combiner C. 

20 in estimating the value of a multi -dimensional 

integral in financial securities trading, a 
multivariate integrand is sampled at points 
corresponding to a low-discrepancy deterministic 
sequence of points in the multivariate unit cube as 

25 defined below. If the multivariate unit cube is also 
the domain of integration, the points of the low- 
discrepancy deterministic sequence can be used as 
sample points directly. In the case of a more general 
domain of integration, sample points correspond to 

30 points of a low-discrepancy deterministic sequence in 
the multivariate unit cube via a suitable 
transformation or mapping. 

When a sufficiently large number of. sampled values 
has been computed, an approximation of the integral is 

35 obtained by suitably combining the computed values, 
e.g., by averaging or weighted averaging. 
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In the d- dimensional unit cube D = [0,l] d , a low- 
discrepancy deterministic s guence z lf z 7 , of 
points in D can be characterized as follows: 

For a point t » [t 1# t d ] in D # define 

5 [0, t) = [0, t,) x ... x [0, t d ), 

where [0, t £ ) denotes a left-closed, right-open 
interval, and denote with Xfo,o(*) the characteristic 

or indicator function of [0, t) . For points z u , 

z B in D, define 
10 R^t; z lf z n ) = (E^, 11 Xp> f o(Zk))/n - tit 2 ...t d# 

and define the discrepancy of z,, . .., z 0 as the L 9 - 

norm of the function R„(.; z,, , z B ) , i,e», 

Hl^t.; z,, . z n ) J„ = supteolRntt; z, f . zj | . 
The sequence z l9 z l9 ... is said to be a low- 
15 discrepancy deterministic sequence provided 

|R.{.; z u . .., z„) \ m m 0((log n) d /n) # 
Low- discrepancy deterministic sequences are 
described in the published literature; see, e.g., 
H. Niederreiter, "Random Number Generation and Quasi - 
20 Monte Carlo Methods", CEMS-NSF . 63, SIAM, Philadelphia, 
1992 and S. Paskov, "Average Case Complexity of 
Multivariate Integration for Smooth Functions", Journal 
of Complexity . Vol. 9 (1993), pp. 291-312. Well-known 
examples of low- discrepancy deterministic sequences are 
25 the so-called Hammersley points, Hal ton points, Sobol 
points, and hyperbolic- cross points. 

For illustration, in the case of Sobol points in a 
single dimension (d=l) , a constructive definition may 
be given as follows: Choose a primitive polynomial 
30 P(x) « x* + ajx"- 1 + ... + a^x + 1 

(whose coefficients a t are either 0 'or l) and define 
so-called direction numbers v j# i > n by the following 
recurrence formula: 

Vi = a^, © a 2 v w 0 . . . © a^v^ © v^ © (v^/2») , 
35 where © denotes a bit-by-bit "exclusive or" operation. 
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The initial numbers v t « m x /2 , . . . , v B «n\/2 a can be 
chosen freely provided m^ is odd and 0 < < 2' for 
i = 1 , 2 , . . . , n. 

Using the direction numbers v { so defined, now 
5 define the one -dimensional Sobol sequence x ( , Xj, . 
by 

= b,v t © bjV, © . . . © b,^, k * 0 
where k « Ejoo^^b^* is the binary representation 
of k. 

10 For higher dimensions (d > 1) , the 

first d primitive polynomials p 1# P 2 , . p d are 
used. If {xfc}^!* denotes the one -dimensional Sobol 
sequence generated by the polynomial P if the d- 
dimensional Sobol points are defined as « (x* 1 , 
15 xf). 

While this definition can be implemented as a 
computer algorithm directly, faster techniques are 
known which produce these points in a "shuffled" or 
permuted sequence. In particular, this applies to the 
20 computer algorithm given in the Appendix. 

For specificity in the following, a preferred 
embodiment of the invention is described as applied to 
a collateralized mortgage obligation known as CMD FN, 
89-23. This has thirty-year maturity and consists of 
25 the following tranches: 

PAC tranches 23-A, 23-B, 23rC, 23-D, 23-E 
supporting tranches 23-6, 23 -H, 23 -J 
residual tranche 23 -R 
accrual tranche 23 -Z 
30 The monthly cash flow is divided and distributed 

according to pre-specif ied rules which cure included in 
a formal prospectus. Some of the basic rules may be 
stated as follows: 

First from the monthly cash flow, the coupon is 
35 paid to the tranches. The remaining amount, called 
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Principal Distribution Amount, is used for repayment of 
the principal. Prior to a fixed date in the future, 
the Principal Amount will be allocated s quentially to 
the tranches 23 -G, 23 -H, 23 -J and 23 -Z. After that 
5 date, the Principal Distribution Amount will be 

allocated sequentially to the tranches 23-A, 23 -B, 23- 
C, 23 -D and 23 -B according to a planned schedule. Any 
excess amount of the Principal Distribution Amount over 
the planned schedule will be allocated sequentially to 
10 the tranches 23-6, 23-H, 23-J and 23-Z. A distribution 
of principal of the tranche 23 -R will be made only 
after all other tranches have been retired. 

In deriving an estimate for the present value of 
the security at the time of issue, the following 
15 notation is used below: 

C - the monthly payment on the underlying 
mortgage pool; 

i k - the projected interest rate in month k, 
k » 1, 2, 360; 

20 «k - the percentage of mortgages prepaying in 

month k; 

a 3«>*+i - the remaining annuity after month k. 
A remaining annuity a* can be expressed as 3 to 
3k -,l + v, + . t . + v 0 w , k - l f 2, . . . , 360, 
25 with v 0 = l/(l+io), where i 0 is the current monthly 
interest rate. Thus, after k months, the remaining 
amount of principal borrowed is C-a^ 

It is assumed that the interest rate i k can be 
expressed as 
30 it - Ko exp<* k )i u , 

where exp(.) denotes exponentiation and where 

%2> — r Ijw ar © independent, normally distributed 
random variables with mean 0 and variance a, *nri Kq is 
a given constant. For the present example, a = 0.0004 
35 is chosen. 
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It is assumed further that w k as a function of i k 
can be comput d as 

w k m Kj + Kj arctan (K9i k + KJ , 
where K u Kj, K 3 , K« are given constants. 
5 Under these assumptions, the cash flow in month k, 

k = 1, 2, — , 360 is 

C (l-w t (£ ,) ) . . . (l-w w (f , $ M ) ) (l-w k (£ t , . . . , £ k ) 

where 

10 VkUt, . f k ) = + K2 arctan (KjKo k i 0 e3qp . .+{ k )+IQ) . 

This cash flow is distributed according to the 
rules of FN, 69-23. Then, the cash flow for each of 
the tranches is multiplied by the discount factor 

. . .v k (£ l# . . . ,f k ) , with 
15 VjC*!,...,^) - 1/(1 + Krfioexp <$|+.. ) , j-l,2,...,360 ff 
to find the present value for month k. Summing of the 
present values for every month gives the present value 
FV X , for each tranche T. 

The expected or estimated value, 
20 E(PV T ) = B(FV T ($„..., 

upon a change of variables is represented by 

E(PV T ) = 1 D PVrty^x,), •.• # y J60 (x^))dx l ...cbC3 60# 

where 

Xi - (2w)-*I^ exp(-t 2 /(2a))dt. 
25 Thus, for each tranche T, a 360-variate integrand 

has to be integrated over the 360 -dimensional unit 
cube, . 

After generating a point 

(X|, .x^, . • • , X^q) 
30 of a low-discrepancy deterministic sequence in the unit 
cube, the point 

(Yi» Yi> Yi6o) 
is computed by finding the value of the inverse normal 
cumulative distribution function at each Xt. Then, for 
35 each tranche T, the function value 
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is confuted. These are the function values used in 
numerical integration. 

Fig. 2 shows results from trial runs for 
5 CMO FN, 89-23 with a preferred method using Sobol 
points generated by the corresponding computer 
algorithm given in the Appendix, as compared with Monte 
Carlo integration. Two Monte Carlo computations were 
carried out, with different "seeds" or starting values 

10 of a congruent ial pseudo - random number generator known 
as RAN2; see W. Press et al.. Numerical Recipes in c. 
Cambridge University Press, 1992. It is apparent that 
the preferred method reaches a steady value more 
rapidly, in this and other trial runs, with typical 

15 requirements of precision and confidence, a speed-up by 
a factor of 3 to 5 was realized as compared with Monte 
Carlo integration. Much greater speed-up can be 
expected when higher precision or/and higher levels of 
confidence are sought. 

20 m a further trial run with CMO FN, 89-23, 

instead of Sobol points, Hal ton points were used as 
generated by the corresponding computer algorithm given 
in the Appendix. It is felt that Sobol points may be 
preferred over Halton points for integrals of high 

25 dimension. However, this preference may not apply in 
the case of lower- dimensional integrals, e.g., with 
dimension up to 5 or so. 

A computation as described may be terminated after 
a predetermined number of function evaluations. 

30 Alternatively, e.g., after every function evaluation or 
after a predetermined incremental number of function 
evaluations, a current approximation may be compared 
with one or several preceding approximations, for 
termination once a suitable condition depending on the 

35 difference betw en approximations is met. Such ter- 
mination criteria may be called automatic. Automatic 
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termination is particularly reliable where a sequence 
of approximations settles down smoothly; see, e.g., 
the curve in Fig. 2 corresponding to Sobol points* 
Advantageously in computing f miction values, a 
5 cluster or network of multiple parallel processors or 
workstations can be used* This may involve a master or 
host processor providing points of a low-discrepancy 
sequence to slave processors and combining function 
values returned by the slave processors into an 

10 approximate value for the integral. Thus, the compu- 
tation can be speeded up in proportion to the number of 
processors used. 

Advantageous further, in combination with 
numerical integration as described above, is the use of 

15 error reduction techniques analogous to variance reduc- 
tion in Monte Carlo integration as described, e.g., by 
M. Kalos et al., Monte Carlo Methods . John Wiley & 
Sons, 1986. This may involve a change of variables 
or/and variation reduction, for example. 
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void sobol (int n) 

i 9 ^** function sobol for generating Sobal points. It returns the 
d : d^s«Mional Sobol point. (The point (0,0, 0) is skipped.) 

The point is implicitly returned through the array x. The function 
sobeeq from Bumsrical Recipes, IB92 was used as a basis, but 
changed significantly to accomodate the parallel distributed approach. 
Ttm function aobol can generate Sobol points skipping an initial 
part of the sequence. The constant MAXDXM, see Numerical Recipes, 
is extended to 360 which required adding more initialising data to 

****** ip (the imitive polynomials), 
£SLil?!f r < ^ grecs i' *** lv (thB initial direction numbers) . The 
polynomial x is used to generate the first coordinate of Sobol points. 

•include "nrutil.h* 
«def ine MAXBXT 30 
ftdef ine MAXDXM 360 



extern int d; /* actual dimension of the points */ 

extern double •x; /* This returns implicitly the n th Sobol point in x */ 
void soboKint n) 
i»t j,k,l; 

unsigned long i,im,ipp; 
static double fac; 

static unsigned long in, ixtMAXDIM+l) , *iu (KAXBXT+1) ; 
static unsigned long mdeg {HRXDXM+l) -(0,M«BIT,1,2,3,3,4,4, 

5,5,5,5,5,5, 

6, 6, €,€,«, 6, 

/♦I**/ 7 ' 7 * 7 ' 7 » 7 ^,7,7 # 7 # 7,7,7,7,7,7.7,7,7, 
' ' . ^ . •f**«»S»« f 8,8,B,e,8,8,B,B,8,8,8, 

/*20*/ in ,« ,„ ,J'?' 9 :?' 9 ' 5 ' 5 ' 9 '5^,9.9,5,9,S,9, 

10 in J2^S r J2^2'JS'^j; aO ' 10 ' 10 ' 1O ' 1O ' 1O '«^O,10,lO,lO t iO,10, 

i2'«'J2 # J2'^^S^2'J2' 10 ' xo ' 18 ' 10 ' 10 - 10 '«,io,io,io io io 

i^i^l^i^l^?^^^^^^^J^^^^ l ' lx ' ll * ll • ll • ll • ll • lI •« 

ii'iz •rr' 11 ' 11 * 11 * 11 ' xx * 11 » 11 »ii»n»ii#ii.ii.ii.ii.ii n n 

11.11.11.11,11,11,11 11 IX, IX, 11, II, 11. 11. II It n it ii ii 

11.11.11.11.11.11.11,11,11,11,11,11 11 ii,il ',11 ii ii 11 11 
".".ii. ii.ii.ii, 11.11.11.11,11,11 ii ii iLiI S'ii ii S' 

".12:i^ l2 '"• la • l^ ' la • la • la • la ' t2 •" : "'"'"'«'«'".«. 

/• The fimt coordinate in , et to MAZBXT, but it i- never u.«d •/ 
static unsigned long iplMWn>XM*ll.{0.»aBIT.0.1.1,2.1,«. 
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/♦dsg 5*/ 2,13.7, 14,11. 4, 

/♦deg «•/ 1,1C,13,22,19,25, 

/♦tog 7*/ l f 32, 4, 8, 7, 56, 14, 28, 19, 50, 21. 42, 31 ,62, 37, 41, 55, 59, 

/♦deg •♦/ 14,56,21,22,38.47,49,50,52.67,70,84,97,103,115,122, 

/♦deg 9*/ 8,13 ,16 ,22 ,25,44 ,47 ,52 ,55 ,59 ,62 ,67 ,74 ,81 ,82 ,87 ,91 ,94, 
103 . 104 , 109 , 122 , 124 , 137 , 138 , 143 , 145 , 152 , 157 ,167,173, 176 , 181 , 182 , 185 , 
191,194,199,218,220,227,229,230,234,236.241,244,253, 

/♦deg 10*/ 4, 13, 19, 22, 50, 55, 64, 69, 98, 107, 115, 121, 127, 134, 140, 145, 1S2, 

158,161,171,181,194,199,203,208,227,242,251,253,265,266,274,283,289,295, 
301,316,319,324,346,352,361,367,382,395,398,400,412,419,422,426,428,433, 
446, 454 , 457, 472, 493 , 505, 508 , 

/♦degll*/ 2,11,21,22,35,49,50,56,61,70,74,79,84,88,103,104,112,115,117,122,134, 
137,146,140,157,158,162,164,168,173,185,186,191,193,199,213,214,220,227, 
236,242,251,256,259,265,266,276,292,304,310,316,319,322,320,334,339,341, 
345,346,362,367,372,375,376,381,385,388,392,409,415,416,421,428,431,434, 
439,446,451,453,457,458,471,475,478,464,493,494,499,502,517,510,524,527, 
555 , 560 , 565 , 569 , 578 , 580 , 587 , 589 , 590 , 601 , 607 , 611 , 614 , 617 , 618 , 625 , 62S , 635 , 
641 , 647 , 654 . 659 , 662 , 672 , 675 , 682 , 684 , 669 , 695 , 696 , 713 , 719, 724 , 733 , 734 , 740 , 
747,749,752,755,762,770,782,784,787,789,793,796,803,805,810,815,824,829, 
830,632,841,847,849,861,871,878,689,892,901,908,920,923,942,949,950,954, 
961 , 968 , 971 , 973 ,979 , 982 , 966, 998 , 1001 , 1010 , 1012 , 

/*degl2*/ 41.52,61,62,76,104,117,131,143,145,157,167,171,176,181,194,217,236, 
239,262,283,286.307}; 

static unsigned long iv (HAXDXH*HAXBIT+11 • { 
0, /* 40 elements in a row */ 

1,1,1,L,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1*1*1.1.1*1»1»1>1*1>1^« 
l,l.l,l,l,l,l,l,l,l,l,l,l,l,l,l,l # l,l,i,l,l,l,l,i,l,l,l,l»l*X«l*lrl*lf 1,1,1,1,1, 
l,l,l.l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,i,l,l,l,l»l#l.l»l»l,l,l,l,l#l»l»l#l# 

l,l,l.l,l,l,l,l,l,l,l,l,l.l.l,l,l,l,l,l,l,l.l,l,l.l,lrl»l,l,l»l»l.l»l»l»l#l-l.l- 

l,l,l,l,l,l,l r l,l,l.l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l,l.X,X.l,l,l,l,l,l,l#l,l,l. 
1,1. l.l,l,l,l,l,l,l,l,l,l,l,l,l t l,l,l,l,l,l,l,l,l,l,l,l,X, 1,1, 1,1, 1,1, 1,1, 1,1.1, 
l,l,l,l,l,l,l,l,l,l,l.l,l,l,l,l.l,l,l,l,l,l,l,l,l,l,l#l#1.3»l»ltl»l#l»l»l»l.l»l» 
1,1,1,1.1,1,1,1,1,1,1,1,1,1,1,1,1.1,1,1,1,1,1,1,1.1,1,1,1,1«1»1»X, 1,1, 1,1, 1,1,1, 
i,l,i,l,i,i,i,i,l,l,l,l,l,l,l,l,i,l.i,i,i,i,i,i,i,i,l,i,i,l,l.l,l#i.l.lrl,l,l,l, /*eod 1 

•/ 

1,3,1,3,1,3.1,3,3.1,3,1.3,1,3,1,1,3,1,3,1,3,1,3,3,1,3,1,3,1,3,1,1,3,1,3,1,3,1,3, 

1,3.1,3,1,3,1,3,3,1,3,1,3,1,3,1,1,3,1,3,1,3,1,3,3,1,3,1,3,1,3,1,1,3,1,3,1,3,1,3, 

1,3,1.3,1,3,1,3,3,1,3,1,3,1,3,1,1,3,1,3,1,3,1,3,3,1,3,1,3,1,3,1.1,3,1,3,1.3,1,3, 

1,3,1,3.1.3,1.3,3,1,3,1,3.1.3,1,1,3,1,3,1,3,1,3,3,1,3,1,3,1,3,1,1,3,1,3,1,3,1,3. 

1,3,1,3,1,3,1,3,3,1,3,1,3,1,3,1,1,3,1,3,1,3,1,3.3,1,3,1,3,1,3,1,1,3,1,3,1,3,1,3, 

1,3,1,3,1,3,1,3,3,1,3,1,3,1,3,1,1,3,1,3,1,3,1,3,3,1,3,1,3,1,3,1,1,3,1,3,1,3,1,3, 

1,3,1,3,1,3,1,3,3,1,3,1,3,1,3,1,1,3,1,3,1,3,1,3,3,1,3,1,3,1,3,1,1,3.1.3,1,3.1,3, 

1,3,1,3,1.3,1,3,3,1,3,1,3,1,3,1,1,3,1,3,1,3,1,3,3,1,3,1,3,1,3,1,1,3,1,3,1,3.1,3, 

1,3,1,3,1,3,1,3,3,1,3,1,3,1,3,1,1,3,1,3,1,3,1,3,3,1.3,1,3,1,3,1,1,3,1,3,1,3,1,3, /♦end 2 

*/ 

1,5,7.7,5,1,3,3,7,5,5,7,7,1,3,3,7,5,1,1,5,3,3,1,7,5,1,3,3,7,5,1,1,5,7,7,5,1,3,3, 

1,5,7,7,5,1,3,3,7,5,5,7,7,1,3,3,7,5,1,1,5,3,3,1,7,5,1,3,3,7,5,1,1,5,7,7,5,1,3,3, 

1, 5, 7, 7, 5,1, 3, 3, 7, S, 5, 7, 7, 1,3, 3, 7, 5, 1,1, 5, 3, 3, 1,7, 5, 1,3, 3, 7, 5, 1,1, 5, 7, 7, 5, 1,3, 3, 

1,5,7,7,5,1,3,3,7,5,5,7,7,1,3,3,7,5,1,1,5,3,3,1,7,5,1,3,3,7,5,1,1,5,7,7,5,1,3,3, 

1,5,7,7,5,1,3,3,7,5,5,7,7,1,3,3,7,5,1,1,5,3,3,1,7,5,1,3,3,7,5,1,1,5,7,7,5,1,3,3, 

1,5,7,7,5,1,3,3,7,5,5,7,7,1,3,3,7,5,1,1,5,3,3,1,7,5,1,3,3,7,5,1,1,5,7,7,5,1,3,3, 

1.5,7,7,5,1,3,3,7,5,5,7,7,1,3,3,7,5,1,1,5,3,3,1,7,5,1,3,3,7,5,1,1,5,7,7,5,1,3,3, 

1,5,7,7,5,1,3,3,7,5,5,7,7,1,3,3,7,5,1,1,5,3,3,1,7,5,1,3,3,7,5,1,1,5,7,7,5,1,3,3, 

1,5,7,7,5,1,3,3.7,5,5,7,7,1,3,3,7,5,1,1,5,3,3,1,7,5,1,3,3,7,5,1,1,5,7,7,5,1,3,3, /♦end 3 

*/ 

1,15,11,5,3,1,7,9,13,11,1,3,7,9,5,13,13,11,3,15,5,3,15,7,9,13,9,1,11,7,5,15,1,15,11,5,3,1,7, 
9, 

1,15,11,5,3.1,7,9,13,11,1,3,7,9,5,13,13,11,3,15,5,3,15,7,9,13,9,1,11,7,5,15,1,15,11,5,3,1,7, 
9, 
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J;».".«.*.».7.t # lJ.U.1.1.7.f.S,XJ.M.XX.J.U.I.a.X».7.».XJ.t.X.XX.7.S.W.X.M.aX.I.3.X.7. 
X.X5.XX.S,3.X.7,9,XJ.XX.X,3.7,9.S.X3.13.XX.3,XS.5.3.XS.7.9.X3.9.X.XX.7.S.X5.X.X5.X1.S.3,X.7. 
l.X5.XX.5.3.X.7.9.i3.XX.1.3.7.9.S.X3.X3.Xl,3.X5,5.3.15.7.9.X3.9.X.11.7.5.X5.X.X5,XX.5.3.X.7. 
X.X5.XX.5.3.X.7.9.X3.XX.1.3,7.9. 5 .X3,X3.XX.3.X5,5.3.X5.7.9.X3.9.X.1X > 7.S.X5.X # X5.XX,5.3.X.7. 
X.X5.XX.5.3.X.7.9.X3.X1.X.3.7.9.5.X3.X3.1X.3.X5.5.3.XS.7.9.X3.9.X.1X.7.5,X5.X.X5.XX. 5 .3.X.7. 
X.X5.XX.5.3.X.7.9.13.1X.1.3.7.9.5.X3.U.X1.3.X5.5.3.X5.7.9.X3.9.X.XX.7.5.X5.X.X5.XX.S.3.X.7. 

^"i"* 7 ^ 31 - 9 ' 3 ' 27 '^ 

^"^"if' 31 ' 9 ' 3 '"' 15 ' 29 ' 21 - 23 ' 19 ' 11 - 25 - 7 ' 13 -"^ 

^"xs 3 . 7 ^?' 31 ' 9 ' 3 - 27 ' 15 ' 29 ' 21 - 23 ' 19 ' 11 - 25 - 7 '"-^ 

^"ii 3 , 7 ^,?- 31 - 9 ' 3 ' 27 - 15 ' 2 ^ 

j:"^",?' 31 ' 9 - 3 ' 27 ' 15 ' 29,21 ' 23 '"'"-"' 7 -"'" ^.».».».31.XX. 5 .23.37.19.2X.S.X. l 7.X 

i:"^"^ 31 ' 9 ' 3 ' 27 -^^ 
^"x"^";?' 31 ' 9 ' 3 ' 27 ' 15 - 29 ' 21 ' 23 ' 19 ' 11 ' 25 ' 7 '"-" 
^"ift 7 *"^ 31 ''' 3 ' 27 '^^ 
^"x&'j"^ 11 /!;^^ 

t":«:;^'?5:i?;^:i2: 35 * 53 ' X9 - 51 '«' 3 ^ 
l:ll:U:Z:Z:%:il^ 

S:ll:'l:ll%:lVd^ 

i:S:£:2:!i\-g:^ 

l^:tl:"^l:l 9 7 :tl:ll^ 

5:S:S:S:S:5^^ 
ls* s <tl:ilt¥&^ 
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1,85, 67. 49,125,25,109, 43. 89, €9,113, 47, 55,97,3 , 37, 83. 103,27, 13,33,115,41,79,17, 29. 119, 7S, 73,1 
05,7,59,(5,21,3,113,61,89,45,107, 

1, 85, 67. 49, 125. 25, 109, 43, 89, 69, 113, 47, 5S, 97, 3, 37, 83, 103, 27, U, 33, 115, 41, 79, 17, 29, 119, 75, 73,^ 
05,7,59, 65,21,3,113, 61, 89, 4S.107, 

1, 85, 67, 49, 125, 25,109, 43, 89, 69, 113, 47, 55, 97, 3, 37, 83,103,27, 13, 33, 115, 41, 79, 17, 29, 119, 75, 73,1 
05,7,59,65,21,3,113,61,89,45,107, 

1,85,67.49,125,25,109,43,89,69,113,47,55,97,3,37,83,103,27,13,33,115,41,79,17,29,119,75,73,1 
05,7,59,65,21,3,113,61.89,45*107, /*end 7 */ 

1,255,79,147,141,89,173,43,9,25,115,97,19,97,197,101,255,29,203,65,195,177,105,17,47,169,109 
,149,15,213,135,253,191,155,175,63,47,7,23,39, 

1,255,79,147,141,89,173,43,9,25,115,97,19,97,197,101,255,29,203.65.195,177,105,17,47,169,109 
,149,15,213,135,253,191,155,175,63,47,7,23,39, 

1,255,79,147,141,89,173,43,9,25,115,97,19,97,197,101,255.29.203,65,195,177,105,17,47,169,109 
,149,15,213,135,253,191,155,175,63,47,7,23,39. 

1,255. 79, 147, 141, 89, 173, 43. 9. 25, 115, 97, 19, 97. 197, 101, 255. 29, 203. 65. 155,177,105, 17, 47, 169, 109 
,149,15,213,135,253,191,155,175,63.47,7,23,39. 

1,255,79,147,141,09,173.43.9.25,115,97,19,97,197,101,255,29,203,65,195,177,105,17,47,169,109 
,149,15,213,135,253,191,155,175,63,47,7,23,39, 

1,255,79,147,141,89,173,43,9,25,115,97,19.97,197,101,255,29,203,65,195.177,105,17,47,169,109 
, 149 , 15 , 213 , 135 ,253 , 191, 155, 175 , 63 . 47 , 7 , 23 , 39 , 

1,255,79,147,141,89,173,43,9,25,115.97,19.97,197,101.255,29,203,65,195,177.105,17,47,169,109 
,149,15,213,135,253,191,155,175,63,47,7,23,39, 

1,255,79,147,141,89,173,43,9,25,115,97,19,97,197,101,255,29,203,65.195.177,105,17,47,169,109 
,149, 15, 213, 13S, 253, 191, 155, 175, 63, 47,7,23, 39, 

1,255,79,147,141,89,173,43,9,25,115,97,19,97,197,101,255,29,203,65,195,177,105,17,47,169,109 
,149,15,213,135,253,191.155,175,63,47,7,23,39, /*«* 8 */ 

1,257,465,439,177.321,181,225,235,103,411,233,59,353,329,463,385,111,475,451,263,19,249,275, 
369,393,167,333,473.469,101,21,451,229,247,297,403,497,61,361, 

1,257,465,439,177,321,181,225,235,103,411,233,59,353,329,463,385,111,475,451,263,19,249,275, 
369,393,167,333,473,469,101.21,451,229,247,297,403,497,61,361, 

1,257,465,439,177.321,181,225,235,103,411,233,59,353,329,463,385,111,475,451,263,19,249,275, 
369,393.167,333,473,469,101,21.451,229,247.297,403,497,61,361, 

1,257,465,433,177.321,181,225,235,103,411,233,59,353,329,463,385,111,475.451,263,19,249,275, 
369,393,167,333.473,469,101,21,451,229,247,297,403,497,61,361, 

h5 5 Zl! 65 ' 09 ' 177 ' 321 ' 181 ' 225 ' 235 ' 103 ^^ 

369,393,167,333,473,469,101,21,451,229,247,297,403,497,61,361, 
h2 5 Zi! 6 V 35 ' 177 ' 321 ' 181 ' 225 ' 235 '"^ 

369,393,167,333.473,469,101,21,451,229,247,297,403,497,61,361, 

h5 5 I ' I 65 * 439 ' X77 ' 321 ' 181 ' 225 • 235 ' 103 • 411 ' 232 ' 59 • 353 * 329 • «« ' «5 , 111 , 47S ,451,263,19,249.275, 
369,393,167,333,473,469,101,21,451,229,247,297,403,497,61,361, 

h2 5 I;i";i 3 ?: r 7 ' 321 • 181 ' 225 

369,393,167,333,473,469,101,21,451,229,247,297,403,497,61,361, 

h2 5 Ii! 6 h! 3 ! ' 1 11 ' 321 ' 18X ' 225 • 235 ' 103 ' 411 * 233 * 59 • 353 ' 329 ' 463 ' 385,111,475,451,263.19,249,275, 

369,393.167,333,473,469,101,21,451,229,247,297,403,497,61,361, 

/*cnd 9 */ 

1 ill X l 721 ' 1013 ' 759 ' 835 ' 949 ' 113 • 929 ' 615 ' XS7 > 39 • 761 ' 169 ' 9 83 ' 557 ' 647,581,505,833,139.147.203.81 
,337,829,989,375,365.131.215,733.451,447,177,57,471,979,197,251, 

1 ;ZI l i! 2X I X0X3 * 759 ' 835 ■ 949 ' 1X3 • «9 , 61S , 1S7 , 39 , 761 , 169 , 983 , 657 , 647, 561 , 505 , 833 , 139 , 147 , 203 , 81 
,337,829,989,375,365,131,215,733,451,447,177,57,471,979,197,251, 

1 ;ZZ 1 ;Z! 1 ii2 1 2 * 759 • 835 • 949 • 1X3 • 929 • 615 • 157 * 39 • 761 • 169 • 983 • 657 • 647 » 581 • 5 °5 . 833 , 139 , 147 , 203 , si 

,337,829,989,375.365,131,215,733,451,447,177,57,471,979,197,251, 

\V*\ll ^i!! 1 !!! 5 !; 835 9 *** ' 112 ' 929 ' 6X5 ' 157 ' 39 • 761 ' 169 ' 983 ' 657 ' 647 ' 5W.305, 633, 139, 147,203,81 
.337,829,989,375,365,131,215,733,451,447,177,57,471,979,197,251, 

X 4ZZ X lZ5 X l52 X ?iZ^? * 83 ^ * 9 * 9 * ^^3, 925 , 615 , 157 , 39 , 761 , 169 ,983 , 657, 647, 581, 505 , 633 , 139 , 147 , 203, 81 
,337,829,989,375,365,131,215,733,451,447,177,57,471,979,197,251, 

^ZZ 1 ;!!^? 1 ^ 5 ! i 835 ' 949 ' 1X3 ' 929 ^ ^ 

,337,829,989,375,365,131,215,733,451,447,177,57,471,979.197,251, 

: 3 2Z: 8 Z?: 9 ! 9 : 379 ' 3€5 ' X31 ' 215 ' 733 ' 45X ' 447 ' X77 ' 57 ' 47l »579. 197.251, 
^ZZ 1 ;? 21 ;" 1 ?;! 5 ?^ 

,337,829,989,375.365,131,215,733,451,447,177,57,471,979.197.251, 
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1,771, 721,1013. 759,835,949,113,929,615,157,39,761,169,983,657,647,581,505,833,139,147,203,81 

,337,829,989,375,365,131,215,733,451,447,177,57,471,979,197,251, 

/•end 10 */ 

1,1285, 823,727,267, 833,471, 1601, 1341, 913, 1725, 2021,1905, 375,893, 1599, 415, 605, 819,975,915, 171 

S, 1223, 1367, 663. 629, S25, 469, 981, 1667, 1587, 1251, 451, 481, 721, 483, 1209, 1457, 415, 1435, 

1, 1285, 823, 727, 267, 833, 471, 1601, 1341,913, 1725, 2021, 1905, 375, 893, 1599,415, 60S, 819, 975, 915, 171 

5, 1223, 1367, 663, 629, 525, 469, 981, 1667, 1587, 12S1, 451, 481, 721, 483, 1209, 1457, 415, 1435, 

1, 1285, 823, 727, 267, 833, 471, 1601, 1341, 913, 1725, 2021, 1905, 37S, 693, 1S99, 415, 605, 819, 975, 915, 171 

5,1223,1367,663,629,525,469,981,1667,1587,1251,451,481,721,483,1209,1457,415,1435, 

1,1285, 823, 727, 267, 833, 471, 1601, 1341, 913, 1725, 2021, 1905, 375, 893, 1599, 41S, 605, 819, 975, 91S, 171 

S , 1223 , 13 67 , 663 , 629 , 525 , 469 , 981 , 1667 , 1587 , 1*51 , 451 , 481 , 721 , 483 , 1209 , 14 57 , 415 , 1435 , 

1 , 12 85 , 823 , 727 , 267 , 833 , 471 , 1601 , 1341 , 913 , 1725 , 2021 , 19 05 , 375 , 893 , 1S99 , 4 IS , 60S , 819 , 9 75 , 91S , 171 

5,1223,1367,663,629,525,469,981,1667,1587,1251,451,481,721,483,1209,1457,415,1435, 

1,1285,823,727,267,833,471,1601,1341,913,1725,2021,1905,375,893,1599,415,605,819,975,915,171 

S, 1223, 1367, 663,629, 525, 469, 981, 1667, 1587, 1251, 451, 481, 721,483, 1209, 1457, 415, 1435, 

1 , 12 85 , 623 , 727 , 267 , 833 , 471 , 1601 , 134 1 , 913 , 1725 , 2021 , 19 05 , 375 , 893 , 1599 ,415,605,819,975, 915 , 171 

5 , 1223 , 1367 , 663 , 629 . 525 , 469 , 981 , 1 667 , 1587 , 1251 , 451 , 481 , 721 , 483 , 1209 , 1457 , 415 , 1435 , 

1,1285,823,727,267,833,471,1601,1341,913,1725,2021,1905,375,893,1599,415,605,819,975,915,171 

5 , 1223 , 1367 , 663 , 629 , 525 , 4 69 , 981 , 1667 , 1587 , 1251 , 451 , 461 , 721 , 4 83 , 1209 , 1457 , 415 , 1435 , 

1, 1285, 823, 727, 267, 833, 471, 1601, 1341,913, 1725, 2021, 1905, 375, 893, 1S99, 415, 605, 819, 975, 915, 171 

5,1223,1367,663,629,525,469,981,1667,1587,1251,451,401,721,483,1209,1457,415,1435, 
/*end li */ 

1,3655, 4091,987, 1839, 4033, 2515, 579, 3863,977, 3463, 2909, 3379, 1349, 3739, 347, 387. 2881, 2821, 1B73, 
1959 , 1929 , 2389 , 3 251 , 1149 , 243 , 3609 , 1131 , 1701 , 143 , 1339 , 3497 , 2499 , 1571 , 983 , 4021 , 1625 , 3217 , 1163 , 
2977, 

1,3855,4091,987,1839,4033,2515,579,3863,977,3463,2909,3379,1349,3739.347,387,2681,2821,1873, 
1959,1929, 2389, 3251, 1149,243,3609, 1131, 1701, 143, 1339, 3497,2499, 1571, 983, 4021, 1625,3217,1163, 

hli* 5 ' 4091 ' 987 ' X839 ' 4033 • 251S * 579 • 3863 > 977 ' 3463 • 2905 • 3379 > 1349,3739,347,387,2081,2821,1873, 
1959 , 1929 , 2389 , 3251 , 1149 , 243 , 3609 , 1131 , 1701 , 143 . 1339 , 3497 , 2499 , 1571 , 983 , 4021 , 1625 , 3217 , 1163 , 
2977, 

hli SS • 409X ' 987 ' 1839 ' 4033 ' 2515 ' 575 • 3863 • 577 • 3463 • 2909 ' 3379 ' 13 « • 3739,347,387,2881,2821,1873, 
2977 ' 1929 * 2389 ' " 51 ' ' 243 ' 3 609 * 1X31 ' X701 ' " 3 ' 1339 ' 3497 * 2499 ' X571 ' 983 • 4021 ' X625 • 3217 ' 1163 • 

h25"li° 91 ' 987 • 1839 • 4033 ' 2515 ' 579 ' 3863 ' 977 • 3463 ' 29 09 • 3379 ' 1349,3739,347.387,2881,2821,1173, 
2977 ' 1929,2389 ' 325X • 1X49 ' 243 ' 3609 * XX3X ' X70X * 143 • 1339 • 3497,2499,1571,983.4021,1625,3217,1163, 

hX 9 !i!2 9 h? 87 * X839 • 4033 ' 25X5 ' 579 • 3863 ' 977 ' 3463 ' 2909 • 3379 • 1349 • 3739 > 347 , 387 , 2881 , 2821 , 1873 , 
2977 ' * 9 3251 ' 1149 ' 243 ' 3609 * X131 9 1701 ' 143 ' X339 ' 3497 ' 2499 ' 1571 # 983 ' 4021 * X 625 ' 32X7 ' XXM ' 
J' 3 ;55,4091,987, 1839 , 4033 , 2515 , 579 , 3863 , 977 , 3463 , 2909 , 3379 , 1349 , 3739 , 347 , 387 , 2881 , 2821, 1873 , 
£$7 ' X929 ' 2389 ' 3251 ' 1X49 • 243 • 36 09 ' 1X31 • 1701 ' 143 ' 1339 ' 3497,2499,1571,983,4021,1625,3217,1163! 

h«2*f ;t2 9 h?! 7 : X839 • 4033 • 25X5 > 579 • 3863 • 977 * 3463 • 2909 * 3379 • 1349 • 3739 ■ 347 . 387 , 2asi , 2821 , 1873 , 

1959. 1929. 2389. 3251, 1149, 243, 3609, 1131, 1701,143, 1339. 3497, 2499, 1571,983,4021! 162513217; 1163 I 

h2Pf:!Ph! 87 ' 1839 • 4033 ' 25X5 ' 579 ' 3863 ' 977 ' 3463 ' 2909 * 3379 • 3739, 347, 3B7,2B81, 2821,1873, 

1959. 1929, 2389, 3251, 1149,243, 3609, 1131, 1701,143, 1339,3497, 2499, 1571,983,4021, 1625,3217ili63, 
2977 /*end 12 */ 

if tn < d) I 

£or(j«12;j<«MAXBlT;j+*> iT(i*j*MRXDXHI-l; /* Initial lie all direction 

nunbers for the first 
coordinate to 1 */ 
for <j«l,k-0;j<«HAXBrr;j^,k>«MRXDHI) iulj] - tivfk] ; 
for (k.l;k<«MAXDIM;k*+> { 

for (ji=l;j<rtrfeg[)cj iu[j)(k] (MAXBZT-j); 

for tj*mdcgtkj+l; j<-MWtBIT;j++) { 
ipp-iptkj ; 
i=iutj-mdegtkl) fkj ; 
i (i » ndeg(kl) ; 
for (l«mdeg(k) -1;1>»1;1--) ( 
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if (ipp & 1) i *m iu(j-UCX), 
ipp >>■ 1; 

iuCj] [k]«i; 



£ac-1.0/(lL « HAXBXT); 
in-O; 

} 

else 

/* Check if the <n-l> -th number waa generated in the previous call 
to sobol. If not, update in and ix */ 

if(ini*n-l) { 
unsigned long gx*y; 

/• Set ix to 0 */ 

for <k-l/k<-IMIH<d,MWa>IM) ;k++) ix(k).0; 
in-n-1; 

grey»in*(in»l) ; /* Find gray code of in */ 
for (j«l?j<«MAXBIT;i+*) { 

if (gray&i) { /* Only digits which are 1 are used •/ 

ia-(j-l|nCAZDZH; 

for Oct ;k<»XMIH(d,KAXDX2f) ix[k] ivtiink); 

^ gray»-l; 

Ltain; /* Calculate the next vector in the sequence */ 

for (j.lf j<-MMCBIT;j++) { /• Find the rightmost sero bit */ 

if (Mia t 1)) break; 

ia »• 1; 

If tj > HAXBIT) nrerror("MAZBZT too small in sobseq*) ; 
in-(j-l) *HAXDIMr 

for (k-lfk<-Hmi(d,MAXDIM) ;k>4-) ( 
ix(k) A . iv(is»k); 
x(k-l).ix(k)*fac; 

ln*+; 



iundef HAXBIT 
tundef HAZOZH 
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void haltonfint n> 

This is the function haltoo for generating Balton points. 

It returns the n-th d-dimensional Hal ton point. The point ie illicitly 

returned through the array x. The last two digits of n-1 in base p[j) 

are kept in ql [j] and q2[j) . When both digits become p[j)-l, the 

radical inverse function is conputed again. That way the accumulation 

of round-off error ie avoided. In practice, there are not any upper bounds 

on the values of d and n. 



extern int d; /* actual distension of the points */ 

extern int *ql,*q2 f /* ql(j) is the last digit of n-1 in base p[j), q2(j] is 

the digit before the last one */ 
extern double *x; /* This will contain the Ralton point •/ 
extern int *p; /• the first d prim numbers •/ 
extern int *p_l; /* first d prime numbers minus 1 •/ 

extern double *incxl, *incr2; /* incrifj] is l/pljl and l/(p[jl*p(j)) */ 

double find.fi (int p, int n) j /* See below •/ 

void haltonCint n) 

double a; 
int j,nn; 

static int ins_n; /* The default value of ins_n is 0 •/ 

/♦ Oieck if the (n-l)-th number was generated in the previous call 
to h al t o n . if not, update qi, q2, and x •/ 

if (ina_nl.n-l) 



ins_asn-i; 
^ for <j«0; j<d? 

ql[J)«insjn%ptj); 
Q2Cjl- Uns_n/p(j] )%ptf ) ; 
j xtj!.findjCi<p|j) ,ias_n) ; 

ins_n4-* ; 

for(j=0;jcd;j++) 

{ if Cqltjj<p_lljl) 

/* It is easy to update when the last digit is less than p[jj -l */ 

j x{j) «x(j] ♦inert (j); 
^ else if (q2(j]<p_l(jj) 

/* This is the case when the laet digit is pfj]-l and the digit 
before the last one is less than pljl-l */ 

qltjl-0; 

^ xCj)*x{j)+incrl(jJ+incr2[j]-l.O; 
else 



SUBSTITUTE SHEET (RULE 26) 



WO 97/07475 



PCT/US95/10363 



< 17 

/* This ia the case when the last digit ia ptj) -1 and the digit 
before the last one is also ptjl -1 */ 

qMjl-0; 
o2 lj).0; 

nn«n/<pljl*plj))f 
a-nn%p(j); 
nn-nn/p[j] ; 

if (m) x(jj-(a+find fitptj] ,nn))*incr2 (j)*incrl(j] ; 
else xCj)-a*incr2{jj**incrl{j) ; 

. '* 

find_fi(lnt p, int n) 

This returns the radical inverse function fi (p,n) at n for the prime p. 

** ~ — / 

double find.fi (int p, int n> 

int p2,nn, 
double s,fi,incr; 

incr«X.0/ps 

p2-p*p; 

nn-n/p2; 

fi-0.0; 
8-incr; 

while (im > 0) 

fi«-« (nn%p)*s; 

UttaDn/p; 

^ »*»incr; 

/* The two largest components of f 1 are added later to avoid possible 

loss of precisian */ 
fi.<(n/p*p| + fi) / p2; 
fi+«in%p>/ (double) p; 
return fi; 
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Claims 

1 1. A computer method in financial securities trading, 

2 for producing an approximate value for an estimated 

3 value of a financial security, comprising: 

4 providing the computer with financial security 

5 data; 

6 producing the approximate value by numerically 

7 integrating a multivariate integrand whose multi- 
B dimensional integral over a domain of integration 
9 represents the estimated value, wherein numerical 

10 integration comprises: 

11 evaluating the integrand at points in the 

12 domain of integration corresponding to points of a 

13 low- discrepancy deterministic sequence, and 

14 combining the integrand values to produce the 

15 approximate value; and 

16 producing the approximate value as an output 

17 for inspection or/and further processing. 

1 2. The computer method of claim 1, wherein combining 

2 the integrand values comprises averaging. 

1 3. The computer method of claim 1, wherein the number 

2 of integrand values is predetermined. 

1 4. The computer method of claim 1, wherein the number 

2 of integrand values is determined automatically. 

1 5. The computer method of claim 1, further comprising 

2 allocating integrand evaluations among a plurality 

3 of processors. 

1 6. The computer method of claim 1, further comprising 

2 application of an error reduction technique. 
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1 7. The computer method of claim 6, wherein error 

2 r duct ion comprises a change of variables. 

1 8. The computer method of claim 6, wherein error 

2 reduction comprises variation reduction. 

1 9. The computer method of claim 1, wherein the low- 

2 discrepancy deterministic sequence comprises Sobol 

3 points. 

1 10. The computer method of claim 1, wherein the low- 

2 discrepancy deterministic sequence comprises Hal ton 

3 points. 

1 11. The computer method of claim 1, wherein the low- 

2 discrepancy deterministic sequence comprises 

3 Haranersley points. 

1 12. The computer method of claim 1, wherein the low- 

2 discrepancy deterministic sequence comprises 

3 hyperbolic -cross points. 

1 13. The computer method of claim 1, wherein financial 

2 securities data comprise derivative instrument data. 

1 14. The computer method of claim 1, further comprising 

2 using the approximate value in offering the security 

3 for sale. 

1 15. The computer method of claim l, further comprising 

2 using the approximate value in deciding whether to 

3 buy, sell or hold the security. 

1 16. A computer system for financial securities trading, 

2 for producing an approximate value for an estimated 

3 value of a financial security, comprising: 
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4 means for providing the computer with financial 

5 security data; 

6 means for producing the approximate value by 

7 numerically integrating a multivariate integrand 

8 whose multidimensional integral over a domain of 

9 integration represents the estimated value, 

10 comprising: 

11 means for evaluating the integrand at points in 

12 the domain of integration corresponding to points of 

13 a low-discrepancy deterministic sequence, and 

14 means for combining the integrand values to 

15 produce the approximate value; and 

16 means for producing the approximate value as an 

17 output for inspection or /and further processing. 

1 17. The computer system of claim 16, wherein the means 

2 for combining the integrand values comprises 

3 averaging means. 

1 18. The computer system of claim 16, wherein the number 

2 of integrand values is predetermined . 

1 19. The computer system of claim 16, further comprising 

2 means for automatically determining the number of 

3 integrand values. 

1 20. The computer system of claim 16, further comprising 

2 means for allocating integrand evaluations among a 

3 plurality of processors. 

1 21. The computer system of claim 16, further comprising 

2 means for applying an error reduction technique. 

1 22. The computer system of claim 21, wherein the means 

2 for applying an error reduction technique comprises 

3 means for a change of variables. 



SUBSTITUTE SHEET (RULE 26) 



WO 97/07475 



PCT/US95/10363 



21 

1 23. Th computer system of claim 21, wherein the means 

2 f r applying an rror reducti n technique compris s 

3 means for variation reduction. 

1 24. The computer system of claim 16, wherein the means 

2 for producing a low- discrepancy deterministic 

3 sequence comprises means for producing Sobol points. 

1 25. The computer system of claim 16, wherein the means 

2 for producing a low- discrepancy deterministic 

3 sequence comprises means for producing Hal ton 

4 points. 

1 26. The computer system of claim 16 , wherein the means 

2 for producing a low- discrepancy deterministic 

3 sequence comprises means for producing Hammersley 

4 points. 

1 27. The computer system of claim 16, wherein the means 

2 for producing a low- discrepancy deterministic 

3 sequence comprises means for producing hyperbolic - 

4 cross points. 
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